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Abstract 

We investigate unconventional superfluidity in a gas of Fermi atoms with an anisotropic p-wave 
Feshbach resonance. Including the p-wave Feshbach resonance as well as the associated three kinds 
of quasi-molecules with finite orbital angular momenta L z = ±1,0, we calculate the transition 
temperature of the superfluid phase. As one passes through the p-wave Feshbach resonance, we 
find the usual BCS-BEC crossover phenomenon. The p-wave BCS state continuously changes 
into the BEC of bound molecules with L = 1. Our calculation includes the effect of fluctuations 
associated with Cooper-pairs and molecules which are not Bose-condensed. 

PACS numbers: 03.75. Ss, 03.75.Mn, 03.70.+k 
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The search for p-wave superfluidity is the next big challenge in a trapped Fermi gas, after 
the discovery of s-wave superfluidityin 40 K and 6 Li|l|, . Recently, p-wave Feshbach 

resonances have been observed 0, 7,0]. The discovery of p- wave superfluidity will be the first 
realization of pseudo-spin triplet superfluidity in a Fermi atomic gas, quite different from the 
recently discovered singlet s-wave superfluidity jl|, 2I 0, U, 0|- Since the pairing interaction 
associated with a Feshbach resonance can be tuned by varying the threshold energy 2v of 
the resonance (see below), one can probe the p-wave BCS-BEC crossover. The superfluidity 
will continuously change from p-wave BCS-type to a BEC-type of bound molecules with a 
finite angular momentum L — 1, as one passes through the Feshbach resonance. 

As a useful first step, we calculate the superfluid phase transition temperature T c over the 
entire p-wave BCS-BEC crossover regime. We explicitly include the p-wave Feshbach reso- 
nance and associated molecules with three values of L z = ±1, 0. To describe the BCS-BEC 
_fl Q, it i S ™ to inch.de the fluctuation, in the three p-wave Coopet- 
channels and their coupling due to the strong pairing interaction associated with the Fes- 
hbach resonance. Our work is a generalization of Ref. We consider both a single- 
component Fermi gas (where a Feshbach resonance occurs in the same hyperfine state) as 
well as a two-component Fermi gas (where a Feshbach resonance occurs between different hy- 



m\ 



perfine states). The Feshbach resonances in both cases have been recently observed| 
We deal with both a narrow and a broad Feshbach resonance. Although we mainly consider 
a uniform gas in this letter, we discuss T c in a trapped gas in the BEC limit. 
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3. 



p-wave superfluidity in trapped Fermi gases was discussed in the BCS regime at T c 
Very recently, the p-wave gap equation for the order parameter of the superfluid phase in 
the c_ re glo n wa, s o,ve d in Re, Q at T = 0. An attracts interaction in the L + 
partial wave channel was considered. An interesting phase transition in a two-dimensional 
Fermi gas has also been predicted!^]. In contrast to these recent papers, our starting 
point explicitly introduces the molecules which form as a result of the Feshbach resonance. 
We thus emphasize the physical nature of the p-wave bound states which form the Bose 
condensate in the crossover region. We also remark that a p-wave pairing mechanism has 
been proposed using the dipole interaction 

We extend the coupled fermion-boson (CFB) model for a s-wave Feshbach resonance 3, 
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la . Il7| to a p-wave one, 
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Here c|, is a creation operator of a Fermi atom with the kinetic energy e p = p 2 /2m. b^j 
describe three kinds of molecular bosons (labelled by j = x, y, z), all with the center of mass 
momentum q, associated with the p-wave Feshbach resonance. The threshold energy 2v in 
the molecular kinetic energy + 2v = q/2M + 2v is independent of j, due to the spherical 
symmetry of the system we are considering. In the last term, g r is the coupling constant of 
a p-wave Feshbach resonance, with pj characterizing the p-wave symmetry 18]. 

The Feshbach resonance term in ((TJ is obtained from a more general Hamiltonian i^F.R. 
j drdr'\g T (r—r')$(r,r')w(r)w(r')+h.c.]. Here \l/(r) = Z) p e 2p ' r c p is a fermion field operator, 
and 

$(r,r')=E eiqR E u nL {f)Y L , Lz {9A)Kn,L,L z (2) 

Q n,L,L z 

describes molecules with center of mass R = (r + r')/2, and relative coordinate f = r — r'. 
^q,n,L,L z is an annihilation operator of a bound molecular state, described by the eigen- 
function u n L{r)YL t L z (0,(j)). The last term in (^J) is obtained when we retain the terms in 
\P(r) to leading order in p, in the L = 1 channel (& q L z = &q,n,L=i,L z ) for a Feshbach res- 
onance state specified by a radial quantum number n. We note that p x oc Yn + ii _i, 
p y oc Y11 — ^1,-1 and p z oc Yi )0 . Thus the molecular operators 6 qj - in (JTJ) are re- 
lated to & q ,L z , with azimuthal angular momentum components L z = ±1,0, as follows 
(6q,a» 6q,v» & q,«) = (^[ & q,i + & q -1 L 73 [ & q,i - &q-i],&q,o)- Equation (JTJ) also includes a non- 
resonant p-wave interaction [/ 18(, which we take to be attractive {—U < 0). 

In the p-wave Feshbach resonance, since two Fermi atoms form one of the three kinds of 
quasi-molecular bosons described by fe q • (j = x, y, z) and this bound state can dissociate 
into two Fermi atoms, we take M = 2m and impose the conservation of the total number of 
Fermi atoms as N = N-p+2 Y^j=x,y,z -^b- Here Np is the number of Fermi atoms and is the 
number of Bose molecules in the j-th component. This constraint can be actually absorbed 
into (JJJ) by considering the grand-canonical Hamiltonian H = H — fiN. The resulting 



Hamiltonian has the same form as (0), where e p and are replaced by £ p = e p — // and 
£q = £q — 2/i, respectively. 

The superfluid phase is characterized by three anisotropic p-wave Cooper-pairs Aj(p) = 
U J2 P > PjP'ji c -p' c p') an d three molecular BEC order parameters 4>j = (6 q=0j ) (J = 
x,y,z). In the equilibrium state, these are related to each other through the identity [lfl] 
Pj^j = 2^-2/^ i(p)- ^he single-partid 6 excitations have the BCS spectrum E p = 



£p + I Z)j A,-(p)| 2 with the composite order parameter, given by Aj-(p) = Aj(p) — y/2g T pj(j)j. 
The angular dependence of A 3 -(p) is proportional to p^. This composite order parameter 
is self-consistently determined by the BCS gap equation, 1 = |?7 e ff J2 P tanh where 
the factor p 2 /3 comes from the angular integration of p|. The effective pairing interaction 
U e g = U + <7j / (2z/ — 2/i) includes the effect of Feshbach resonance ljj. In the weak-coupling 
or BCS regime, \l ~ e-p (where ep is the Fermi energy). 

The analogous p-wave CFB model for a two-component Fermi gas (=f, |) is described by 



~ U H P " P /c P +q/2T C -p+q/2i C -p'+q/2| C p'+q/2T 
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In a mean field pairing approximation, we again obtain the same single-particle excitations 
E p and the gap equation as those in the single-component case. The Cooper-pair order 
parameter is A 3 -(p) = U J2 P ' PjP'j( c - P > i c p'l) ■ 

We now present the p-wave strong-coupling theory at T c for the single- component (spin 
polarized) model defined in (HI). The discussion is easily extended to the two-component 

nn 

case. The equation for T c is obtained by employing the Thouless criterion |9|, [10|, the tem- 
perature when the particle-particle scattering matrix first develops a pole at to = q = 0. 
In the t-matrix approximation, the p-wave scattering matrix has the form T^p, p', q, uo) = 
PiTij((i,u))pj, which is shown diagrammatically in Fig. 1(a). In this figure, the first and the 
second lines, respectively, describe the effects of non-resonant interaction U and the p-wave 
Feshbach resonance, that give f (q, w) = {IY,-} = —[1 — U e g(q, w)f[(q, uj)]~ 1 U e g(q, u) = 
x,y, z). Here, £/ e fj(q, u) = U — g^Do(q,u) is an atom-atom interaction including dynamical 
effects described by the bare molecular Bose propagator Dq 1 (q, to) = uj + iS — [£q + 2v\ . The 
correlation functions II = {11^} are obtained from the analytic continuation iv n — > + i5 
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of the two-particle thermal Green's function, 

n„(q, w„) - \ £ P , ft ^%f^%^, (4) 

P p ?P+q/2 + <5p-q/2 - W n 

where f(e) is the Fermi distribution function. The diagonal components ITjj (i = x,y,z) 
describe superfluid fluctuations in the i-th Cooper- channel, while the off-diagonal compo- 
nents give the coupling of fluctuations in different channels. Noting that 1X^(0,0) = in 
our approximation, the Thouless criterion gives the equation for T c as 



This has the same form as the mean-field gap equation at T = T c , with Aj(p) — > 0. 
However, the chemical potential u in (jSJ) can be quite different from the Fermi energy e F in 
the crossover regimej^, Q and one needs an additional equation to determine /i. 



The chemical potential fi is determined from the equation for the number of atoms N, 
which is calculated from the thermodynamic potential Q using the formula N = — ^jjj], 
[nj. Figure 1(b) shows the fluctuation correction to Q, where the diagrams on the left 
and right describe fluctuations in the p-wave Cooper-channels and the Feshbach resonance, 
respectively. Summing up these diagrams, we obtain 

N = N ? ~ 4 E ^-tr\\ogD(q,iu n ) 

^-tr[\og[l-UIL(q,ip n )] 
= N F + 2N B + 2N C , (6) 

where the trace is taken over the L = 1 space (j = x, y, z). N F = J2 P /(£p) is the number of 
free Fermi atoms. is the number of Feshbach molecules, given as the poles of the renor- 
malized (matrix) molecular Bose Green's function .D _1 (q, ii/ n ) = iv n — + 1v) — S(q). 
The molecular self-energy E(q, iu n ) = —q?TL/(l — UIl) describes the fluctuation effects in 



the p-wave Cooper-channels. As in Refs. 



3, liy], N c can be interpreted as the contribu- 



tion of preformed p-wave Cooper-pairs as well as particle-particle scattering states. We 
note that superfluid fluctuations in the three p-wave Cooper- channels are strongly coupled 
to one another in (0) through (i ^ j). Equations (0) and (JBJ) are the basic coupled 
equations describing T c of a uniform p-wave superfluid over the entire BCS-BEC crossover. 
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The same equations are obtained in the two- component case in (J3J), with Np replaced by 
Np = 2^ p /((p), reflecting the two Fermi hyperfine states. 

Figure 2 shows the self-consistent numerical solutions of (j5J) and (JBj). In this figure, we 
have introduced a p-wave scattering length a„ for the renormalized interaction U3c which 

n n 

occurs in the gap equation when written in a cutoff-independent way 21 . This is defined 11, 
12j by -4ir(3a 3 p )/m = = U eS /{l - E[o,u c ] |r)> where u c is a high-energy cutoff. 



The increase of (kpa p )~ 3 corresponds to a decrease of bare threshold energy 2v. Since 
the chemical potential also decreases to approach v [see the inset in panel 2(b)], the bare 
interaction U c g = U + g^/(2u — 2/x) becomes stronger for larger (kpa p )~ 3 . In the BCS regime 
[(kpa p )~ 3 < — 1], T c agrees well with the standard weak-coupling BCS theory ['BCS' in Fig. 
2(a)]. On the other hand, in the crossover regime [—1 < (kpa p )~ 3 < 0], the deviation of 
T c from the weak-coupling result is large. The chemical potential fi also begins to strongly 
deviate from the Fermi energy Ep, as shown in Fig. 2(b). Figure 3 shows that the gas 
continuously changes from a gas of Fermi atoms (dominated by Np) into a Bose gas of 
bound states (dominated by Am = Ab + Nc). In the BEC regime [(kpa p )~ 3 < 0], free Fermi 
atoms are almost absent, and T c approaches a constant value. Its precise value depending 
on whether one is dealing with a single-component gas (ft) or a two-component gas (TJ.)- 
This difference is due to different Fermi energies in the two cases (see TABLE I). The peak 
in T c in Fig. 2(a) would be absent if the coupling to the bound states [h], |3] was properly 
included in the self-energies of the Fermi atoms. 

In the extreme BEC limit, where all the atoms have formed Feshbach molecules 
(A^f, Ac = 0), the gas can be regarded as a non-interacting Bose gas mixture with 
three kinds of Feshbach molecules, with L z = ±1, 0. In this case, rewriting Eq. © as 
2/i = 2u - g?Tlii{0, 0)/[l - UTLu(0, 0)], we find 2/x -> 2v, because IL^O, 0) = in this BEC 
limit. This result is consistent with the inset in Fig. 2(b). Since 2/x is the chemical poten- 
tial of the molecular Bose gas and 2v is the threshold energy of molecular excitations, the 
condition 2/x = 2v is that required for BEC in a non-interacting Bose gas. That is to say, 
T c in this extreme case is simply determined by A^ = 3^ q nB(£q), where ub{s) is the Bose 
distribution function. The factor 3 comes from the presence of three kinds of molecules, 
which is characteristic of p-wave superfluidity. Because of this factor, T c in the p-wave case 
is lower than the s-wave case, as shown in TABLE I. TABLE I also shows T c in the BEC 
limit in a trapped gas, evaluated within the LDA. These values for T c in a trapped gas 
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seems accessible in current experiments. The crossover behavior of T c shown in Fig. 2(a) is 
a general result valid for any type of p-wave superfluidity. 

Figures El and El indicate that the crossover behavior of T c , /i, and the number of Fermi 
atoms Np and that of Bose molecules Am show quasi-universal behavior when plotted as 
a function of (kpa s )~ 3 , irrespective of whether the Feshbach resonance is narrow (g r < 
£p) or broad (g T > ep). On the other hand, the character of the bound state bosons is 
different between the two. In a narrow Feshbach resonance, the Feshbach molecules (N-q) 
are dominant in the crossover regime, while Cooper-pairs (iV c ) are dominant in a broad 
Feshbach resonance (see Fig. 3). However, Feshbach molecules always dominate in the 
extreme BEC limit. 

We note that the phase diagrams in trapped Fermi gases which experiments 

measure involve passing through the resonance in an adiabatic (constant entropy) manner. 
In the case of a s-wave Feshbach resonance, Ref. j^j] has discussed how one can calculate 
such phase diagrams using a simple ideal gas model. This could be extended to the p-wave 
Feshbach resonance case. 

To summarize, we have discussed the BCS-BEC crossover in the presence of a jo-wave 
Feshbach resonance. Generalizing earlier work on the s-wave BCS-BEC crossover [10], we 
have included fluctuation effects in the three p-wave Cooper-channels, as well as the three 
kinds of Feshbach molecules with L z = ±1,0. Observation of the molecular condensate in 
the BEC regime |l|, Q. [1 U, 0] would be a first step in the study of p-wave superfluidity. 

I would like to thank Prof. A. Griffin for discussions and also critical reading of this 
manuscript. This work was financially supported by a Grant-in-Aid for Scientific research 
from Ministry of Education of Japan and funds from NSERC of Canada. 
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FIG. 1: (a) Particle-particle scattering matrix in the ^-matrix approximation in terms of the 
non-resonant interaction U (first line) and the p-wave Feshbach resonance g T (second line). Go 
and Dq are the bare single-particle Fermi and Bose Green's function, respectively, (b) Corrections 
to the thermodynamic potential originating from fluctuations in the p-wave Cooper-channels (left 
diagram) and the Feshbach resonance (right diagram). 

FIG. 2: (a) T c in the p-wave BCS-BEC crossover, jf: single-component Fermi gas, with U = 
NppU = 0.4ep- Tl : two-component Fermi gas, with U = 0.8£f- Results for a narrow Feshbach 



resonance {g r = y Np^g T = 0.6£f) and a broad Feshbach resonance (g r = 5ef) are shown. 'BCS' 
shows a weak-coupling result in the two-component case with fi being fixed at the value at v = 
2.5£f- (b) Chemical potential fi(T c ) in a single component Fermi gas. The solid and dashed lines 
show the results for a narrow and a broad Feshbach resonance, respectively. The inset shows /u(T c ) 
as a function of the threshold energy 2v for a narrow Feshbach resonance. 



FIG. 3: Numbers for various kinds of particles at T c in a single-component Fermi gas. (a) narrow 
Feshbach resonance, and (b) broad Feshbach resonance. Nm = Afe + ^Vcj where iVe describes 
Feshbach molecules and iVc gives the contribution from Cooper-pairs (stable and unstable). 
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TABLE I: T c in the BEC limit. (S) and (T) show the single- and two-component Fermi gas, 
respectively. T c in a uniform gas is given by T c = j- — , 12/3 (where ((z) is the zeta- function), 
with a = 1 (s-wave), a = 6 (S), and a = 3 (T). In a harmonic trap, T c = ^ a ^j 1/3 is evaluated 
using the LDA. Tp (= £f) is obtained from = i]J2 £ < £F 1 5 where 77 = 1 (77 = 2) for the single 
(two) component case. 

symmetry uniform gas [Tp] trapped gas [Tp] 

s-wave 0.218 0.518 

p-wave (S) 0.066 0.285 

p-wave (T) 0.105 0.359 
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